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PREFACE  AND  SUMMARY 

In  view  of  the  paucity  of  experimental  data  on  the  nature  of 
the  electromagnetic  signals  from  nuclear  explosions  and  their 
effects  on  military  systems,  it  has  been  necessary  to  place  consid¬ 
erable  emphasis  on  theoretical  prediction.  Such  theoretical  predic¬ 
tion  on  the  nature  of  the  signals  requires  knowledge  of  the  explosion- 
induced  air  conductivity.  In  the  past,  there  has  been  substantial 
uncertainty  about  how  to  calculate  this  conductivity,  arising  from 
the  uncertainty  in  the  spectrum  of  secondary  electrons  as  they  slow 
down.  In  the  present  report,  this  spectrum  is  estimated,  and  the 
electronic  conductivity  per  electron  is  shown  to  be  equal  to  that  for 
a  mean  electron  energy  corresponding  to  thermal  and  electrical 
equilibrium.  For  sea-level  nuclear  explosions,  this  mean  energy 
is  determined  by  the  equilibrium  electric  field  of  about  1  esu  and  is 


about  0.  16  ev. 
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AIE  CONDUCTIVITY  PRODUCED  BY 
NUCLEAR  EXPLOSIONS 


In  order  to  determine  the  electromagnetic  field  generated  by 
nuclear  explosions  in  the  atmosphere,  it  is  essential  to  specify  the 
explosion— induced  conductivity  of  the  air  in  the  neighborhood  of  the 
burst  point.  It  is  known  that  this  conductivity  depends  upon  the 
spectrum  of  the  ionized  electrons  through  the  dependence  of  the  con¬ 
ductivity  on  the  energy— dependent  electron— atom  collision  frequency. 
In  the  present  report,  the  spectrum  of  the  ionized  electrons  is 
determined  in  terms  of  their  source  function.  With  this  spectrum, 
the  electron  collision  frequency,  and  hence  the  electronic  conductivity 
per  ionized  electron,  is  estimated. 

ELECTRON  CONDUCTIVITY 

The  total  electronic  contribution  to  the  explosion— induced  air 
conductivity  is 

cMt)  =  J  dE  cre(E,  t),  (1) 


where 


a  (E,  t)  =  — 
e  '  m 


n(E,  t) 
v(E)  ' 


(2) 


v(E)  is  the  collision  frequency  of  the  slow  secondary  electrons  with 
air  molecules  and  is  approximately1 

v(E)  ~  3  x  109  p  E  sec  1,  ^ 

1D.  R.  Bates,  ed.  Atomic  and  Molecular  Processes,  New 
York,  1962.  Chap.  10. 
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where  p  is  the  air  pressure  in  mm  of  mercury  and  E  is  the 
electron  energy  in  e.v.  n(E,  t)dE  is  the  density  of  electrons  with 
energies  between  E  and  E  +  dE. 

An  accurate  evaluation  of  n(E,  t)  can  be  made  using.  Age 
Theory — namely,  for  each  collision  of  an  electron  with  air,  the 
electron  is  assumed,  to  undergo  its  mean  energy  loss.  If  /3(E)  is 
the  electron  attachment  rate  at  energy  E,  y(E)  the  recombination 
rate,  and  \(E)  the  fractional  energy  loss  per  electron  collision  with 
air,  then  Age  Theory  gives 


=-pn-7N.n+-|(\Env)+  S(E,  t),  (4) 

where  S(E,  t)  dE  is  the  production  rate  of  electrons  with  energy 
E  to  E  +  dE  per  unit  volume  and  time  and  N.  is  the  density  of 
positive  ions . 

N\  is  related  to  the  electron  density  through  the  equation 


9N. 

1 

9t 


=  -  a  N. 


o  l  L  i 


N .  -  /  dE  n(E 


dE  y(E)  n(E,  t)  + 


dE  S(E,  t) 
(5) 


where  a  is  the  neutralization  rate  for  positive  and  negative  ions, 
o 

Charge  neutrality  implies  that  the  negative  ion  density  is 
N.  -  J dE  n(E,  t). 

A  straightforward  solution  of  Eq.  (4)  gives  the  electron 
density  in  terms  of  the  source  strength  and  the  density  of  positive 


ions: 


— 3— 


n(E,  t) 


1 

X.  E  v 


s(e\  t 


dE"  \ 
ETrXTT"vTr  ) 


(6) 


.  exp 


dE" 


|3(E")  +  y(E")  N.  (t 


r  e" _ dE _ \ 


Accordingly,  the  electron  conductivity  is 


a  (t)  =  — 
e  m 


dE 


dE 


o  EX v 


exp 


1  S^E t  - 


,E' 


dE" 


'E 


E "x"v" 


) 


,.E’ 

-  /  dE" 


/3(E")  +  y(E")  N.  (t-  j 


,  E  " 


(7) 


dE 111 


E 


E"'X"'  v" 


J- 


E 


EX  v 


and  the  total  electron  density  is 
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_E' 

E 


dE" 
E'V'v1' 


(10) 


exp 


E' 

-  /  dE" 

JE 


/  rE  dE1"  \  ^  —  1 
/3(E")+y(E")N.(t-  / 

JE  * 


E"  X"  v' 


Since,  in  general,  N.  depends  on  ng  (by  Eq.  (5)),  this  is 
only  a  formal  expression  for  .  In  many  important  cases, 

however,  the  recombination  rate  yN\  can  be  neglected  compared 
to  the  attachment  rate  /3,  so  that  <^>  depends  only  on  the 
source  function  and  can  therefore  be  evaluated  explicitly. 
APPROXIMATIONS  TO  <t)/ne(t) 

Evaluation  of  cr^t)  requires  specifying  S(E,  t).  Since, 
however,  S(E,  t)  depends  upon  the  design  of  the  particular  explosion 
device,  we  cannot  evaluate  o^t)  in  complete  generality.  We  can, 
however,  obtain  an  accurate  approximation  to  ce(t)/ne(t),  the  elec¬ 
tronic  conductivity  per  electron,  which  is  valid  for  any  specified 

S(E,  t),  provided  electron  recombination  can  be  neglected  compared 
2 

to  attachment. 

For  this  purpose,  we  shall  assume  an  approximate  form  for 
S(E),  t),  neglect  recombination,  and  evaluate  a^OO/n^t)  within 
these  approximations.  Then  we  will  indicate  the  conditions  under 
which  recombination  may  be  neglected  and  show  that  our  result  for 

2 

In  the  Appendix  the  effect  of  retaining  electron  recombination 
is  discussed. 
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is  independent  of  the  assumed  form  for  S(E,  t). 

We  shall  assume  that  the  secondary  electrons  are  all  pro¬ 
duced  at  a  single  energy  E^;  then 

S(E,  t)  =  n(t)  6(E  —  E][),  (11) 

where  n(t)  is  the  total  number  of  electrons  produced  per  unit 
volume  and  time. 

Equation  (10)  then  becomes: 


,E 


1  dE  1  . 
EX.v 


-  h(t 

V  \ 


'l  dE1  \ 
,'v'J 


E'k' 


_E 1 


/  ( - I* 

j3(E')+r(E')N  (t-/  ^irii -ii  )' 


exp 


_E, 


dE' 


;E 


E'k'v' 


(12) 


,.E 


1 o  EX.v 


1  dE  .  /.  rE  1  dE*  \ 

n  V1  JE  EVv') 


exp 


.  E , 


/3(E')+y{E')N.(t-/ 


.E1 


dE" 


E  V 


v,r  )-i 


-1 


dE1 
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EVv1 


Neglecting  yN.  compared  to  0,  this  simplifies  to 
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1  dE 


EXi 


3  1  .  fEl  dE'  \  T  fE  1 

“  7  n  WE  iw)  exp  r  JE 


dE’ 


P(E') 

E'X'v' 


.E 


1  dE 


EX 


E  .  /  fE  1  dE1  x 
v  nV  JE  E'X'v'  ) 


exp 


..E 


;E 


'dE'-^ 


E'X'v' 


(13) 


Replacing  E  by 


.E. 


[~1  P(E') 

JB  6E  ® 


(14) 


as  the  integration  variable, 


dx 


1  1  .  u  fX  dx'  \ 

s  7  n1-/ o  — ) 


r 


.00 

/  dx 
o 


1  [*dx'\-x 


— X 

e 


(15) 


If  the  electron  production  rate  is  slow  compared  to  the 
attachment  rate  /3,  then  for  the  important  values  of  x  in  the 
integration, 


(16) 


dx  e 


o 


,  -x  1 
dx  e  p 


1  dE  1_ 
EX  v  v 


1  dE  ■ 
EX  v  6 


H E1) 

E'X'v' 


E, 


dE' 


E 


flE') 

E'X'v' 


(17) 
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Thus  for  slow  production  rates,  y  is  independent  of  the  source 
rate  and  depends  only  on  the  initial  electron  energy. 

If,  on  the  other  hand,  the  electron  production  rate  is  fast, 
or  comparable  to  the  attachment  rate,  the  approximation  of  Eq.  (16) 
cannot  be  made.  For  this  case,  it  is  convenient  to  assume  that  the 
rate  is  exponentially  increasing;  that  is, 


at. 


n(t)  ~  e  ; 


a  >  S  . 


Then 


f  ,  1  i  -x 

/  dx„  -e  exp<  —  a 

>o  13  v  It, 


f  ,  1  -X 

/  dxj-e  exp 


(18) 


(19) 


Basic  Physical  Data 

The  initial  energy  of  the  secondary  electrons  is  about  10  e.v., 
since  on  the  average  33  e.v.  are  expended  per  ion  pair,  and  about 
one-half  of  this  energy  goes  into  ionization.  The  electrons  are 
slowed  down  by  collisions  with  air  molecules,  losing  some  energy 
by  elastic  collisions  (on  the  average  0.  37  x  10  of  their  energy 
per  collision),  but  most  by  inelastic  collisions  which  excite  vibrations 
(at  the  higher  energies)  and  rotations  of  the  air  molecules.  After 
enough  collisions  the  energy  will  be  reduced  to  thermal  energies 
~  1/25  e.v.  (at  kT  ~  ^77  e.v.).  At  thermal  energies,  the 


—8— 


collision  frequency  is  ~  10  sec 

Data  on  the  average  energy  loss  per  collision  of  slow 

electrons  with  air  is  obtained  mainly  by  swarm  measurements  and 

3 

has  uncertainties  of  interpretation.  Massey  and  Burhop  present 

two  sets  of  data  for  X  as  a  function  of  E  ,  the  mean  energy  of 

the  electron  swarm,  which  disagree  considerably  due  to  differences 

in  experimental  determination  of  the  relation  of  E  and  F/p.  (F  is 

the  electric  field  causing  the  drift  of  the  swarm.)  However,  for 

mean  electron  energies  above  0.2  e.v. ,  both  sets  of  data  agree 
—3 

that  X  >  1.3  x  10  .  Moreover,  recent,  more  reliable  data  quoted 

by  Huxley  and  Crompton  for  very  slow  electrons,  E  <  0.1  e.v. , 

-4 

gives  X  >  4.x  10  .  For  our  numerical  estimates,  we  shall  assume 

*V 

X  =  I.  3  x  10  3  for  E  >  0.  2  e.v. ,  and  X  =  4  x  10-"4  for  E  <  0.  2 

e.v.  Since  in  both  cases  we  take  the  minimum  values  for  X  ,  the 

4 

effect  is  definitely  to  underestimate  the  slowing  down;  that  is,  to 
overestimate  the  effective 

Electron  attachment  to  neutral  oxygen  forming  0  involves 

'  o 

a  three-body  process  for  which  at  sea  level  the  rate  is  about  10 
-1  5 

sec  .  At  thermal  energy  (1/25  e.v.  ) 

Massey,  H.  S.W.  and  E.  H.  S.  Burhop,  Electronic  and  Ionic 
Impact  Phenomena,  Oxford,  1952,  p;  279. 

4 

This  underestimate  makes  even  stronger  our  conclusion  that 
the  effective  conductivity  per  electron  corresponds  to  a  near-equilibrium 
spectrum. 

5 

D.R.  Bates,  ed.  Atomic  and  Molecular  Processes,  New  York 
1962,  Chap.  6.  ’ 
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8  8 

the  rate  is  0.75x10  ;  it  rises  to  ~  1.5x10  at  0.1  e.v., 

8 

then  falls  slowly  to  ~  0.4  x  10  at  1  e.v.  For  numerical  estimates, 

8  •— 1 

we  shall  take  /3  =  constant  «  10  sec  ,  at  sea  level;,  or 

g  2  _ 1 

)3  10  (pi Pq)  sec  at  other  densities. 

The  recombination  rate  of  electrons  and  positive  ions,  y, 

g  3  2 

is  about  10  cm  sec  for  the  dissociative  process 

+ 

e  +  0g  “*■0  +  0. 

In  addition,  at  high  pressures,  there  is  a  three-body  process  in— 

0  0  ““l 

volved.  At  sea  level,  this  has  a  rate  of  the  order  of  10  cm  sec  , 
and,  of  course,  becomes  less  important  at  high  altitudes.  Thus, 

g 

the  term  y]SL  can  be  neglected  compared  to  j3  if 

yN.  «  /3 

or 

10"6  N.  «  108  (p/pQ)2;  N.  <<  IQ14  (p/pQ>2. 

The  positive  ion  density  for  a  typical  nuclear  explosion  can 
be  estimated  by  computing  the  y— ray  energy  deposition,  assuming 
33  e.v.  per  ion  pair  produced.  If  is  the  total  prompt  y— ray 
energy,  and  y  is  the  total  yield  in  kilotons,  the  ion  density  at 
distance  r  in  meters  is 

®See  Appendix . 
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N.  ~  .  6  x  10 17  -X.JL  e-r/300  po/p 

l  2  p  ’ 

r  o 

where  it  is  assumed  that  E^  =  3  x  10~^  Y.  For  N.  <<  10 14  (p/p  )2  , 

Y  -r/300  p  /p  .  in-3  , 

— ^  e  o  <<  1.  6  x  10  p/p  . 

^  o 

r 


For  example,  at  sea  level,  for  a  1  KT  explosion,  the  in¬ 
equality  is  satisfied  at  distances  r>  25  meters.  For  a  1  megaton 
explosion,  the  distance  is  about  400  meters.  At  higher  altitudes, 
for  kiloton  yields  the  distance  increases  as  {pj p)  ;  for  megaton 
yields  at  very  high  altitudes  Qpl  PQ  <  10  3^),r  >  800[Y^,p  pQ/p]1^2 

7 

meters.  Thus,  if  we  consider  distances  greater  than  this  value  of 
r,  we  can  neglect  vN.  compared  to  3,  and  use  the  approximate 


formulae  for  (  — 

v 


NUMERICAL  ESTIMATE  OF  cr  (t)/n_(t) 

_ _ _ e  B 

It  is  a  straightforward  matter  to  substitute  the  energy  de¬ 
pendent  functions  X,  3  ,  and  v  into  the  approximations,  Eqs.  (17) 
and  (19),  and  evaluate  a  (t)  numerically.  However,  it  is  sufficiently 
accurate  to  approximate  these  functions  by 
3  (E)  =  3  =  108  sec  1 

12  —1  —1 

v(E)  =  aE,  a  =  2.  3  x  10  sec  e.v. 

X  (E)  =  X.  a  1.  3  x  10~3,  E  >  E, 
l  b 

=  X2  ~  4  x  10^ ,  E  >  Eb 
-  — ^ 

If  E  >  3  x  10  Y,  the  critical  values  of  r  are  increased. 
For  example,  'With  Ey  =  3  x  10— 2  y,  which  is  a  reasonable  upper  bound 
of  Ev,  r~  200  meters  for  Y  =  1  kiloton  and  r~  1100  meters  for  Y  = 
1  megaton  at  sea  level. 
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where  E^  £  0.2  e.v.  These  approximations  enable  Eqs.  (17)  and 
(19)  to  be  evaluated  analytically. 

For  the  initial  energy  we  take  E^  =  10  e.v.  Then 


fEl  j3  /3 
x  =  /  dE  -  -z- 


JE 


EX  v  aX 


JL. 

aX  \  E  E^ 


,  E  >  E, 


E<  v 

2  b 


where 


(20) 


*.J_  (-1- 

d  aX,  \E; 

1  b 


-TT) 


(21) 


Then 


1 

aE 


aE  +  /3 


i +(x-V  f 

b 


x<  *b 
x>  xb. 


(22) 


Thus  Eq.  (19)  becomes 


\  V 


1  \  _  a  + 

v/=  ~  j3 


r*b  ^(1+f>,  1  /n 

l0  dx*  (lar/r) 


d 


(23) 
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a  +P  [  1 


X1  ^  1  ^  \,a  (' 

T+lE;L+e  1 


E, 


i) 


■i 


(. 


i 


a+/3  j3  a  +fi  VaE 


1 

aE" 


Substituting  numerical  values  for  the  parameters,  and  setting 

g  _ ^ 

a  %  10  sec  ,  we  find  from  Eq.  (19) 


^  ~  0.  6  x  10  * *sec. 


or 


~  1.  7  x  10^  sec  which  corresponds  to  an  average 


electron  energy  of  about  .08  e.v.  This  result  is  essentially  in¬ 
dependent  of  the  initial  energy  E  .  A  somewhat  larger  value  for 

8 

Xj  gives  an  average  energy  of  E^,  independent  of  the  other 

parameters.  For  our  case  this  is  about  0.1  e.v. 

Equation  (17)  leads  to  the  same  result,  with  a  set  equal  to 
zero.  Neglecting  unimportant  terms,  we  find 


/J_\  * 

\  V  /  ~  0 


xrx2  —  aEb 
- 6 


(24) 


or  expanding, 


T/  «  “T  +  l1-^)  5-5xl° 


(■ 

&> 


-12 


sec. 


(25) 


*v  1.  8  x  10^  sec  \ 


which  again  corresponds  to  an  energy  of  ~  0. 1  e.v. 
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The  calculations  above  have  assumed  that  the  electrons  in 
colliding  with  the  air  tend  to  lose  all  their  energy.  Actually,  of 
course,  they  only  tend  to  slow  down  to  ~  0.  04  e.v. ,  the  ambient 
thermal  energy.  Since  the  effective  average  energy  is  about  twice 
as  large  as  this,  no  significant  correction  is  necessary. 

If  there  is  an  electric  field  present,  however,  the  electrons 
at  equilibrium  would  have  a  drift  velocity  which  will  increase  their 
energy  above  thermal.  If  this  new  equilibrium  energy  is  greater 
than  the  0. 1  e.v. ,  which  we  found  above,  then  it  will  determine 
the  air  conductivity. 
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APPENDIX 

It  is  obvious  from  Eq.  (7)  that  if  recombination  to  positive 
ions  is  not  neglected,  then  cr^t)  is  less  than  we  have  estimated. 

This  lessening  results  from  two  effects.  First,  the  electron  density- 
decays  more  rapidly  due  to  the  additional  loss  mechanism.  Second, 
and  perhaps  more  important,  the  lessening  is  greatest  at  the  low 
(near  thermal)  energies  which  tends  to  raise  the  mean  electron 
energy  and  lower 

Since  the  dependence  of  Y  (E)  on  E  is  unknown,  it  is  not 
possible  quantitatively  to  determine  the  change  in  <7  (t)  or 
resulting  from  recombination  effects.  Nonetheless,  it  is  possible 
to  indicate  more  precisely  than  done  in  the  text  the  conditions  under 
which  the  neglect  of  recombination  is  unimportant. 

If  recombination  is  to  be  unimportant,  then,  as  before, 

yn.  <  <  e . 

i 

In  the  text  we  assumed  that  this  inequality  must  hold  for  N.  equal 
to  the  total  number  of  positive  ions  which  are  formed.  Actually  N.  , 
which  obeys  Eq.  (5),  attains  this  value  only  if  S  (t)  is  a  6  — function, 
and  then  only  initially.  A  more  reasonable  estimate  of  N.  can  be 
obtained  from  Eq.  (5)  if  it  is  assumed  that  o^=  Y  (E)  ,  a  constant. 
Experimentally,  a  ~  y  (E  ~  0)  and,  at  higher  energies,  Y  (E) 
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decreases  somewhat.  Under  the  approximate  assumption,  a  y  (E), 
we  find 


SN. 
_ x 

St 


=  -  +  S  (t)  .. 


This  equation  can  be  solved  in  two  limits.  First,  if  S  (t)  =  SQ  6  (t)  , 


N. 

1 

where  S 


o 


if 


S 

o 

1  +a  S  t  ’ 

o  o 

is  the  total  number  of  positive  ions  produced,  and,  second, 

*  0, 


o 


For  the  impulsive  source, 


y  N.  ~  a  N.  <<  P 
l  ox 

implies  that  (c^  ~  y) 


yS 


1  +  yS  t 
o 


<<  P  • 


This  inequality  is  satisfied  for  all  times  in  the  region 
text  and  in  all  regions  about  the  burst  point  for  times 


indicated  in  the 
t  >>  P_1  • 


However,  since  the  production  of  gammas  in  a  nuclear 
explosion  extends  over  times  T  substantially  longer  than  1/P  (P 
evaluated  at  sea  level),  it  is  a  better  approximation  at  sea  level  to 
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assume  that 


Since  T  ~  several  x  10  and  y 


N. 


/102S  o 
(several) 


(at  sea  level), 


Or 

yN.  <<  p 

gives 

14 

S  <<  several  x  10 
o 

at  sea  level.  This  condition  is  somewhat  less  restrictive  than  in 
the  text  and  decreases  the  radius  at  sea  level  outside  of  which 
yN.  «  P  by  a  factor  of  2  or  so. 

Furthermore,  since  N.  decays  as 


Is  / TY 

O 

1 +T,/ 

's  /Ty  (t  -  T) 

O 

after  the  gamma— ray  source  cuts  off  (  t  >  T),  y  N,  <  <  P  is  satisfied 
everwhere  after  times 

t  —  T  >>  P_ 1  . 


